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Efficient Sampling of Band-limited Signals 
from Sine Wave Crossings 

J. Selva 



Abstract — This correspondence presents an efficient method 
for reconstructing a band-limited signal in the discrete domain 
from its crossings with a sine wave. The method makes it 
possible to design A/D converters that only deliver the crossing 
timings, which are then used to interpolate the input signal 
at arbitrary instants. Potentially, it may allow for reductions 
in power consumption and complexity in these converters. The 
reconstruction in the discrete domain is based on a recently- 
proposed modification of the Lagrange interpolator, which is 
readily implementable with linear complexity and efficiently, 
given that it re-uses known schemes for variable fractional-delay 
(VFD) filters. As a spin-off, the method allows one to perform 
spectral analysis from sine wave crossings with the complexity of 
the FFT. Finally, the results in the correspondence are validated 
in several numerical examples. 



I. Introduction 

The analog-to-discrete (A/D) conversion is the first step 
for the discrete-time processing of continuous signals. This 
conversion is fundamentally based on the Sampling Theorem, 
which states that a band-limited signal can be recovered from 
its regularly-spaced samples taken at least at twice the Nyquist 
rate. However, some authors early noticed that this recovery 
is also possible from the signal's zeros, or from its crossings 
with another signal like a sine wave, [l]-[4]. This is due to 
the fact that a band-limited signal is an entire function of 
exponential type, for which there is a factorization in terms of 
its roots akin to that of conventional polynomials, (Hadamard's 
factorization theorem [5, chapter 2]). A consequence of this 
is that it would be possible, in principle, to design A/D 
converters in which the sample quantization is substituted 
by a zero crossing detector and an accurate timing, [6], 
This new procedure would eliminate the need to quantize 
any signal samples, so decreasing the complexity and power 
consumption of A/D converters, provided there is an accurate 
timing available. Besides, it would mainly re-use existing 
technologies, given that zero crossing detection is implicit in 
many existing systems. This last point can be readily seen in 
the current trend in A/D converter design, in which the sample 
amplitudes are turned into zero crossings, which can then be 
accurately detected with low-power consumption, [7], [8]. 

The main obstacle for this alternative procedure is how the 
signal should be reconstructed or processed in the discrete 
domain, since Hadamard's factorization theorem does not 
directly lead to efficient implementations, due to its slow 
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convergence rate. Here, the usual approach in the litera- 
ture consists in approximating the signal in a finite interval 
using a trigonometric polynomial, [6], [9], [10]. But then 
the interpolation error decreases only as 0(l/iV), while the 
complexity per interpolated value is 0(N), where N is the 
number of crossings inside the interval. So, in this approach 
it is necessary to employ a large N to ensure an acceptable 
accuracy, with the associated high complexity. 

The purpose of this correspondence is to present a method 
for overcoming this obstacle, that makes it possible to re- 
construct the bandlimited signal from its sine wave crossings 
efficiently. The method is based on viewing the reconstruction 
as a problem of interpolation from nonuniform samples, to 
which the efficient technique in [11] is applied. Relative to 
the state of the art, it has several advantages: 

• The complexity is reduced significantly. If in the approach 
in [6], [9], [10], a complexity 0(N) (per interpolated 
value) gives an interpolation error 0(l/iV), with the 
proposed method a complexity 0(N) gives an error 
0(e~ 7Ti - 1 ~ BT ' >N ), where B is the signal's two-sided band- 
width and T is the average crossing separation. In practice 
this means that "any" accuracy can be achieved with a 
small 7Y. 

> The method is based on the evaluation of a fixed smooth 
function and on the Lagrange interpolator. Besides, it 
can be evaluated with cost O(N) per interpolated value, 
and can be implemented by re-using efficient designs for 
variable fractional-delay (VFD) filters, [11, Sec. IV]. 

• As a spin-off, the method permits one to perform spec- 
tral analysis from sine wave crossings with complexity 
0(N log N), while the usual method has complexity 
0(iV 2 ), [6, Sec. IV]. 

The correspondence has been organized as follows. The next 
section reviews the state of the art and presents the problem 
formulation. In it, it is shown that the reconstruction from sine 
wave crossings can be turned into an interpolation problem 
from nonuniform samples. Then, this last problem is addressed 
in Sec. [Ell] where the solution adopted is that in the recent 
reference [11], based on the Lagrange interpolator. The main 
result is a simple interpolation formula for reconstructing the 
signal from its sine wave crossings. Then, Sec. [IV] addresses 
the problem of analyzing the spectrum from sine wave cross- 
ings in the light of the formula in Sec. [ill] It turns out that 
this formula makes it possible to reduce the complexity from 
the usual 0(iV 2 ) order to order O(NlogN). Finally, Sec. |V] 
validates the results in the correspondence through a numerical 
example. 

II. State of the art and problem formulation 

The usual representation for a real finite-energy signal s(t) 
with spectrum inside [-B/2, B/2] is given by the Sampling 
Theorem. If the samples s(n/B) are known, this theorem 
states that s(t) can be perfectly reconstructed using the series 

oo 

s(n/B + u)= ^2 s((n-p)/B)smc(p + Bu), (1) 
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where n is an integer and u is any time shift following 
— 1/(25) < u < 1/(25). This series is the basis of most 
processing algorithms for continuous signals in the discrete 
domain. Several authors [1], [2], [4] soon noticed that s(t) 
can alternatively be viewed as a polynomial of infinite degree, 
which can be described in terms of its roots. The reason why 
is that a band-limited signal can be regarded as an analytic 
function over the whole complex plane, if the t variable is 
allowed to take complex values. Besides, this kind of function 
is bounded on the real axis (real t), and its maximum growth 
rate is that of e^ 5 '*' along the imaginary axis. It can be shown 
that this kind of signal admits the representation 

s(t) = K lim TT' (1-t/rp), (2) 

R— >oo 

\t p \<R 

where K is a real constant, and the t p are complex roots which 
appear in conjugate pairs, r p ^ 0, and may not be distinct, 
[12, Theorem VI]. The prime Q means that the factor should 
be replaced with t if t p = Lj. Eq. (0 is the explicit root 
factorization of s(t) as infinite-degree polynomial, and from 
it it is obvious that the roots t p determine the signal except 
for the scale factor K. The root density in any circle \t\ < r 
follows the same rule as the sample density in (Q]), i.e, the 
circle \t\ < t contains either 2r5 samples in (0 or 2tB 
roots in asymptotically. 

The factorization in Eq. (0 suggested the possibility of 
performing the usual processing of continuous signals in the 
discrete domain using the zeros of the signal instead of 
its samples. In such approach, the analog-to-discrete (A/D) 
conversion would consist in acquiring the signal's zeros, and 
the discrete-to-analog (D/A) conversion would be based on the 
evaluation of a formula like 0. However, it was soon realized 
that there were two main problems, [4], [6]. The first was how 
the roots should be located efficiently in the A/D conversion, 
since they may have a non-null imaginary part. And the second 
was how the infinite product in should be approximated so 
as to implement the D/A conversion, since in practice only a 
finite sequence of roots t p is known, and converges slowly. 

A solution for the first problem was readily found [4, Sec 
V], and consisted in subtracting a sinusoidal A s sin(7r5i) to 
s(t), where 

A s = sup |s(t)|. (3) 

t real 

This simple procedure actually solved the problem of locating 
the roots, because the subtraction of this sine wave "moves" all 
roots to the real axis, due to a theorem of Duffin and Schaeffer, 
[13]. Specifically, s(i) — A s sin(7r5i) can only have zeros on 
the real axis, and each of them can be viewed as a zero of 
Asxn(nBt), which has been shifted by at most 1/(25). So, 
in notation, all zeros of s(t) — A s sin(7r5i) have the form 

'The factorization in J5) is due to Titmarsch [12, Theorem VI], and is 
a refinement of Hadamard's factorization theorem [5, chapter 2] for s(t) 
bounded on the real axis. Note that there is a bug in Eq. [4] of [11]. This last 
equation is Hadamard's factorization and its correct form is 

oo 

m = Ae-^ H' (1 - t/t«)e*/*". 



n/B + 5 n , where n is an integer and \S n \ < 1/(25), and two 
consecutive zeros may overlap only at instants n/B + 1/(25). 
Using this description, the factorization in Eq. (0 for s(t) — 
A s sin(7r5t) is 

oo 

s(t) - As sin^ Bt)=K' (t -t) 11 (4) 

fc=i ~ k k 

where K' is a real constant and 

t n = n/B + S n . (5) 

From Eq. @, the basic design for the desired A/D converter 
was clear. First an oscillator would be used to generate the 
wave A s sm(ivt/B), which would then be subtracted from 
s(t). Afterward, the zero crossings would be detected using 
a gate, and the converter output would be the sequence of 
shifts S n in Eq. (0. 

As to the second problem, the usual solution to date 
consists in approximating the signal using a trigonometric 
polynomial, [6]. In short, if the finite sequence of roots 
t n ,t n +i, . . . ,t n +M-i is known, then s(t) — A s sin(7r5t) is 
interpolated using a trigonometric polynomial of order M 
which is zero at t n + m , < m < M. However, this solution 
is not satisfactory since its accuracy is poor even for a large 
number of roots. This complexity issue is the main obstacle 
for achieving efficient implementations. 

The purpose of this correspondence is to provide an efficient 
solution to this second problem. The key point is to realize that 
approximating s(t) from the roots t n of s(t) — A s sin(7r5i) is 
the same as approximating s(t) from its value at these instants, 
since s(t n ) — A s sin(7r5i„). So, this is actually a problem of 
interpolation from nonuniform samples. 

III. Proposed solution 

Reference [11] presents an efficient and accurate inter- 
polator for band-limited signals from nonuniform samples, 
applicable to signals of the same type as s(i), (bandwidth 5, 
supremum amplitude A s ). Consider one such signal z(t) and 
assume it is known at a set of instants r p of the form 

T p = P T + i lp , -P< P <P, (6) 

where the period T and the shifts r) p follow BT < 1 and 
rjp < T/2 (strict inequalities). The approach in [11] consists 
in applying the Lagrange interpolator to the product z(t)j(t), 
where j(t) is a fixed function, and then solving for the value 
of z(t). The interpolator's formula is 

where 

p 

l(*)= n t~Tp, (8) 

p=-p 

and j(t) is given in Ap. Uand only depends on 5, T and P. 
Usually, Eq. (0 is applied for \t\ < T/2, though wider ranges 
are allowed. In [11], it was shown that the accuracy of (0 
increases exponentially with P. So, in practice a small P is 
enough to obtain high accuracy. 
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To assess the accuracy of (0, assume BT = 0.7 and 
\z(t)\ < 1. Then, following the analysis in [11, Sec. Ill], the 
interpolation error of (Q for the function 7(i) specified in Ap. 
U is well fitted by 



e (dB) w 4.12106 + 66.6044 5 - 9.35838 8 2 

- 8.30873P + 3.13419 SP - 0.125803 S 2 P, 



(9) 



where S is a bound on the deviations of the instants r p from a 
uniform grid with spacing T and <$ < T/2. So, if <5 = T/4 the 
value P = 10 gives an interpolation error below e = —55 dB, 
and P = 16 gives an error below e = — 100 dB. Any practical 
accuracy can be obtained by slightly increasing P for fixed 8. 
(For a detailed analysis, see the previous reference.) 

Coming back to the problem in the previous section, the 
sine wave crossings are equivalent to nonuniform samples of a 
form similar to that in ||7). Specifically, the crossing instants in 
© have a maximum deviation from the uniform grid, |<5„| < 
1/(25), and the sample values are given by the sine wave, 
since 



s(t n ) = A s sm(TtBt n ) = A s (-l) n sm(irB6 n ). 



(10) 



However Eq. (Q is not applicable to s(t) since the regular 
grid in (© matches exactly the Nyquist rate and it may 
be \S n \ = 1/(21?), while in (0 there is some sampling 
inefficiency, BT < 1, and it is always \r] p \ < T/2. These latter 
conditions can be easily imposed on the sampling scheme by 
slightly increasing the amplitude and frequency of the sine 
wave. So if A s sin(7rPi) in is replaced by another sine 
wave Asxn.{nt/T), with amplitude A > A s and semi-period 
fulfilling T > and BT < 1, the theorem of Duffin and 
Schaeffer in the previous section ensures that the only zeros 
of the signal 

s(i) - Asm{irb/T) (11) 



are simple and occur at positions 

t n = nT + <L 



(12) 



for integer n. Besides, the shifts 6 n now follow the strict 
inequality |<5„ < S < T/2, where 



5 = (T/ir)axcwx(A s /A), 



(13) 



since there can be no zero crossings whenever 
A\ sin(vrf/T)| > A s . 

To check these conditions numerically, consider the BPSK 
signal in Fig. Q] This signal was generated by modulating a 
raised cosine pulse of roll-off 0.2 and bandwidth B = 0.7/T 
with a sequence of random amplitudes ±1. Then the signal's 
peak amplitude was scaled to 1. Fig. [2] shows the zone 
marked with a rectangle in Fig. [T] together with the sine wave 
Asm(irt/T). The conditions in (fT2l and (Tl3T > mean that the 
crossings with the sine wave can only take place inside the 
shaded rectangles, and there is exactly one in each of them, 
as can be seen in this example. 

Let us derive the final interpolation formula. Assume s(i) 
must be interpolated at an arbitrary t from the crossings 
with Asin(7rf/T), (A > A s , BT < 1). Any t can be 
written uniquely in the form t = nT + u with integer n and 




100 



Fig. 1. Random BPSK signal with bandwidth B = 0.7/T and peak 
amplitude 1. The symbols used to generate this signal were random ±1 values, 
and the modulating pulse was a raised cosine with roll-off factor 0.2. The 
signal's peak amplitude was scaled to one. 
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Fig. 2. Short piece of the BPSK signal in Fig. [TJ overlapped with the sine 
wave A sm(nt/T), with A = 1.1, The stems indicate the positions of the 
crossings with the sine wave. These can only take place inside the rectangles, 
and there is exactly one in each of them. 



-T/2 <u< T/2. Besides the signal z(u) = s(nT + u) with 
time variable u is of the same type as s(t), since a time shift 
affects neither the bandwidth nor the supremum amplitude. So 
(0 is valid on z(u) with instants t p = pT + 8 n + p and values 
z(pT + S n+p ) = A(-l) n+ P sin(7r<J n+p /r), i.e, it is 

1 

s(nT + u) 
p 



7(u) 



E 



A(-l) n +P sin(7T(W P /Tb(pT + S n+p )L n (u) 



p=-p 
_ A(-l) n 



L n(P T + Sn+p)(u -pT- 5 n+p ) 



(14) 



E 



(-l) p sm(Tr6 n+p /T)-f(pT + S n+p )L n (u) 



L 'n(P T + $n+p)(u -pT - 8 n+p ) 
Ln(w) = || I' i - u n+p . 



p=-P 

where 

p 

n u~pT 

p=-p 

Eq. (TT~4b is the final interpolation formula. 



(15) 
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(b) P = 3 

3. Sample BPSK signal in [2] and interpolated signal using I14t for 

: 2, 3. 



For the example in Fig. [2] Fig. |3(a)| shows the signal and its 
interpolated version using (TBI) for P = 2. The discontinuities 
in the interpolated signal take place at t = nT + T/2, integer 
n, because the set of crossings used is different for each n. 
Notice however that just five crossings (P — 2) give a good 
accuracy. Fig. |3(b)| shows the same comparison but for P = 3, 
(seven crossings). The difference between both signals is much 
smaller than in Fig. |3(a)| For P > 3 the error becomes too 
small to be represented this way. For an error analysis see Sec. 

m 

The formula in Eq. ( TBI yields discrete-time processing 
methods for delivering samples of s(i) with any spacing, 
simply by assigning proper values to n and u. The simplest 
case is for spacing T, simply by setting u = 0. For a generic 
grid of instants n\T\ with T\ > and integer n±, the grid 
sample s(niTi) is obtained from ( TBI ), simply by setting n 
and u equal to the modulo-T decomposition of ri\T\ in (TBI) , 
i.e, 



\n\T\jT + 1/2J and u = n{T x - nT. 



(16) 



As to the efficient implementation of Eq. ( TBI ), it was shown 
in [11, Sec. IV] that it can be evaluated with cost just O(N). 
See also the numerical examples in [11, Sec. V]. 
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Fig. 4. Error in the interpolation of s(t) from its sine wave crossings versus 
index P. The error norm is the supremum sup t6 j \s{t) — s(t)| where s(t) 
is the interpolated signal. 



IV. Spectral estimation from sine wave crossings 

The formula in ( TBI makes it possible to interpolate the 
input signal at any instant from its crossings with the sine 
wave, and the cost of this operation is just 0(— loge), where 
e is a bound on the interpolation error. This is because the 
error of ( TBI decreases exponentially with trend c^ 7 'i 1 ^ BT ) p . 
So, if e is set below the working numerical error, Eq. ( TBI ) 
allows one to obtain one sample of s(t) with a small and 
constant computational cost. Therefore, the cost of computing 
TV samples in a regular grid with arbitrary spacing T\ is 
0(N). Once these samples are available, the situation is 
the usual one in which the signal's spectrum is estimated 
from regularly spaced samples, and any of the well-known 
techniques in spectral analysis becomes applicable, [14]. Since 
these techniques are based on the FFT whose complexity 
is O(NlogN), it is clear that the total complexity is also 
O(NlogN). A numerical example is presented in the next 
section. 

V. Numerical examples 

A. Sampling of a BPSK signal 

To validate the results in a specific example, a BPSK signal 
s(t) was generated with the following parameters, 

• Modulating pulse: raised cosine with roll-off 0.2. 

• Random amplitudes equal to ±1. 

. Total two-sided bandwidth B = 0.7/T. 

. Time interval I = [0, (N - l)T] with N = 1024. 
Then several numerical experiments were conducted. 

The first experiment consisted in interpolating s(t) in / from 
its sine wave crossings. The result is shown in Fig. @] for A = 
1.1 and A = 16, where the error norm is the maximum over 
/, that is, if s(t) is the interpolated signal, then the ordinate 
in this figure is 

Bup|s(*)-s(*)|. (17) 
tei 

Notice that this error decreases exponentially with P. Besides, 
the values of S for A = 1.1 and A = 16 are 0.36T and 0.02T, 
respectively, but the error is roughly the same in both cases. 
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Fig. 5. RMS error versus the SNR in the interpolation of the samples s(raT) 
from the sine wave crossings of the noisy realization s(t) + w(i). 



Fig. 7. Amplitude spectrum of sequence s(nT) + w(nT) for nT in /. The 
spectrum's maximum value has been normalized to one. 
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Fig. 6. Maximum difference between the samples s(nT) + w(nT) and their 
interpolated values from sine wave crossings versus the SNR. The error norm 
if that defined in fT9l 



So, the fact that the sampling instants may differ from the grid 
nT (integer n) has a minimal effect on the performance. 

In the second experiment, a white noise process w(t) of 
bandwidth B was added to s(t). Then, s(f)+w(£) was sampled 
at instants nT (integer n) in /, and these samples were also 
interpolated from the sine wave crossings of s(t) + w(t) for 
A = 3 and P = 4 and 9. Fig. [5] shows the resulting root-mean- 
square (RMS) error. The crosses (+) indicate the deviation of 
samples s(nT) + w(nT). The other two curves are the RMS 
errors for P — 4 and P = 9, given by 



' £ |s(nT)-sx(nT)| 2V 



./V 



nTel 



(18) 



where si(i) is the value interpolated from the sine wave 
crossings of the noisy signal s(t) + w(i). The curve for either 
value of P overlaps the sample deviation up to an SNR 
threshold which is fixed by the specific value of P. So, below 
this threshold, the performance is the same if either the signal 
is directly sampled, or if it is interpolated from its sine wave 
crossings. The threshold can be fixed to an SNR as large as 
desired, simply by slightly increasing P, due to the exponential 
dependence of the interpolation error on P. 



B. Spectral analysis 

As to the spectral analysis, it is worth comparing the 
conventional procedure from uniform samples, with the one 
proposed in this correspondence from sine wave crossings. 
In the conventional procedure, the samples s(nT) + w(nT) 
would be delivered by an A/D converter, and then any of 
the existing spectral analysis methods would be applied to 
these data, [14]. And in the proposed procedure, the A/D 
converter would deliver the sine wave crossing timings S n , 
then the uniform samples Si (nT) would be computed using 
Eq. ( fl4l >. and finally the spectral analysis would be the same 
as in the conventional procedure, i.e, it would be performed 
on the samples Si (nT) instead of s(nT) + w(nT). The fact 
is that the result of both procedures would be the same up 
to the numerical accuracy in use. This can be readily seen in 
Fig. [6] in which the error measure is the maximum difference 
between s(nT) + w(raT) and Si (nT), 



sup \s(nT) 

nTel 



■ w(nT) — si(nT)|. 



(19) 



This coincidence is due to the fact that the interpolator in [141 is 
also reconstructing the noise realization w(£), since it is also a 
signal with bandwidth B. Fig.Qshows the amplitude spectrum 
of the sequence s(nT)+w(nT), where the maximum has been 
normalized to dB. If this spectrum were computed from the 
sine wave crossings, the it would differ from that if Fig. |7]by 
the amplitude given in Fig. [8] 

VI. Conclusions 

A method has been presented that makes it possible to 
recover a band-limited signal from its crossings with a sine 
wave. It allows one to design A/D converters which only 
deliver the timing of the sine wave crossings, so allowing for 
a smaller complexity and power consumption in the converter. 
The method is based on viewing the problem as one of 
interpolation from nonuniform samples, to which a recent 
efficient technique is applied. This technique is based on the 
Lagrange interpolator and allows for efficient implementations 
based on current designs of VFD filters. As a spin-off, the 
method permits one to perform spectral analysis from the sine 
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Normalized frequency fT 

Fig. 8. Amplitude of the difference between the spectrum of s(nT)+w(nT) 
and that of si (nT). 

wave crossings with the complexity of the FFT. The method 
has been validated in several numerical examples. 
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Appendix I 

Weight function for the Lagrange interpolator 

Following [11], it is first necessary to define a band-limited 
window function w(<) that approximately selects a finite time 
range. This function is the inverse Fourier transform of the 
Kaiser-Bessel window, 

_ (20) 

smc(j B W T W ) 

where 

B w = 1/T — B and T w = PT. (21) 



Note that in (120b the argument of the sine functions may be 
pure imaginary. In this case, the sine can be evaluated from 
the hyperbolic sine since, for real a, it is 

. , sin(77ra) e~ 7va — e* a sinh(7ra) 

sine (ja) = ) J = = y -^. (22) 

J 7 ™ (2j)(j7ra) ttci 

The weight function ^(t) in Sec. [Ill] is then given by 

_ {-If w(t)L (t) 

lit] = {p\f sinfrt/jy (23) 

where L Q (t) is the Lagrange kernel for the instants pT, \p\ < 
P, 

p 

L„(t) ee Yl t-pT. 

P=-P 
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